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Method of analyzing a data set comprising a tubular structure 



21.12.2000 



The pres»t invention relates to a method of analyzing an object data set 



comprising points in a multi-dimensional space, in which dataset a tubular structure occurs* 
said method comprising the following steps: 
a. Choosing a starting position in or near the tubular structure; 
5 b. Deriving a cutting plane through the tubular structure at the starting position* 

c. Determining a number of points forming part of the surface of the tubular structure in the 

vicinity of the starting position; and 
d Calculating a gradient lo the surface for each of said points. 

10 The intemationdl patent application EP00/09S0S (PHN 17.678) of flie same 

applicant (not yet published) relates to a method of the type motioned above. This method 
generally rdiates to the analysis of a tubular structure in a muiti*dimensional space. 
According to the method described in the international patent application EP00A)950S a self- 
adjusting probe is defined for analysis of the object data set The self-^justmg probe 

15 comprises a sphere and a plane fiirou^ the center of the sphere. The sphere should be 

positioned such that the tubular structure intersects the sphere, at least partially. The plane 
should be oriented orthogonal to the tubular structure. When oriented correctly the self- 
adjusting probe enables semi-automatic shape extraction of tube^Iike geometry* 



20 examined. The object data set notably relates to the density distribution in the object to be 
examined; in that case the data values are the local d^ity values of (a part of) the object to 
be examined. The data values may also relate, for example, to the distribution of the 
temperature or the magnetization in the object The multi-dimensional space is usually the 
three-dimensional space. The data values then relate to a volume distribution of die relevant 

2S property, for example, the density distribution in a volume of the object to be examined. The 
multi-dimensional space may also be two-dimensional. In that case the data values relate to a 
distribution of the relevant property in a plane ibrou^ the object, for example the density 
distribution in a cross-section through the object. 



Such an object data set represents one or more prop^es of the object to be 
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The object data set can be acquired in a variety of ways. The object data set 
notably relates to a padent to be examined. Such an object data set can be acquired by means 
of various techniques such as 3D X-ray rotational angiography, computed tomography, 
magnetic resonance imaging or magnetic resonance angiography. 
5 The method described in the international patent application EPO0/O950S is 

particularly suitable for analyzing the structure of a blood vessel. Several physical 
characteristics of the blood vessel can be determined, such as the diameter thereof, which are 
crucial for diagnosis and a save treatment of for instance a stenosis or an aneurysm. 

Currently the calculation of the above mentioned important step of orienting 
10 the cutting plane is based on a selection of surface points lying inside the sphere at a small 
distance ftom the plane. This is described in the article "Semi-Automatic Shape Extraction 
from Tube-like Geometry" by J. Biuijns. published in B. Girod» G. Greiner, H. Niemami, H.- 
P. Seidel (Eds.). Vision Modelling and Visualization 2000, Proceedings* November 22-24, 
2000, Saarbruecksm. In the example described in the article the surface of a blood vessel is 
15 represented by means of triangles (a known method which will be explained later). However, 
the known selection of surface points is ntha axbitxaiy and may unintradedly include points 
associated with neighboring vessels. This directly leads to inaccuracies in the detenninacion 
of the physical characteristics of the blood vessel under analysis. 

20 The invention aims at providing a method of the type mentioned above which 

allows for a more accurate calculation of the physical characteastics of the tubular structure 
under analysis. 

The method according to the invention thereto further comprises the 

25 charactmzing steps of: 

e) For each point determining a vector from the centre of the wbular structure to said point; 

f) DetOTnining the angle between said vector and the gradient ai said point; 

g) Adding said point to a selection of points if said angle is equal to or smaUer than a 
predetermined ceiling value; 

30 h) Using said selection of points to calculate an orientation for the cutting plane such that 
the direction thereof is as parallel as possible to the lon^tudinal axis of the tubular 
structure at the starting position; 
i) Repeating steps a througjh h f or a new starting position along the tubular structure, if 
necessary. 
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IQ a preferred embodiment the method further comprises the steps of: Defining 
a sphere which is at least partially intersected by the tubular structure; and performing steps e 
through g only for surface points lying inside the sphere. When the radius of said sphere is 
chosen wisely, e.g. slightly larger than the escpected radius of the tubular structure undw 
S analysis, points within said sphere most likely belong to said tubular structure whereas points 
outside the sphere probably don't. Fuithermore the reduction of the number of points taken 
into account during the calculation also increases the calculation speed. 

Jn a further preferred embodiment the method steps e thxough g are performed 
only for points lying at a predetermined maximum distance from the cutting plane. This 
10 allows for a further advantageous reduction of the number of points involved in the 
calculation thus further increasing the calculation speed. 



The inventLon also concerns a computer program to cany out the method 
according to the invention. 

15 

The invention wiU be f urthtf explained by means of the attached dmwing, in 

which: 

Figure 1 shows a cross section thrcu^ two adjacent tubular structures under 
analysis by the meUiod according to the invention. 

20 In figure 1 a cross section is shown through two adjacent tubular structures, 

vessels 1 and 2. A sphere 3 is defined around vessel 1. The center of Ae sphere 3 and vessel 
1 coincide and is indicated as c. Sphere 3 also partly intersects vessel 2. 

According to &e preferred embodiment the method according to the invention 
comprises the following steps. Hrst a starting position is chosen in or near the tubular 

25 structure and a cutting plane is derived through the tubular structure at the starting position. 
These steps are also in more detail described in the international patent application 
£P00/09505t which entire docimient is incorporated herein by reference. 

Next a number of points forming part of the surface of the tubular structure are 
defined in the vicinity of the starting position. In order to illustrate this step three points (also 

30 r^erred to as verdces) on the surfaces of the vessels are indicated. Point 4 is a surface point 
on vessel 1. Points 5 and 6 are surface points on vessel 2. In the relevant art several mediods 
are known to define die surface of a vessel in an object data set. When the object data set 
comprises a volumetric representation of the objea by means of voxels the surface of the 
vessel is defined by a set of boundary vessel voxels. This voxel set can be determined by first 
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separating the vessel voxels from the tissue voxels by means of known techniques, such as a 
'region growing algorithm'. Next it is detemiined which of the vessel voxels lie adjacent to 
tissue voxels. These are referred to as boundary vessel voxels. Preferably this detenninaiion 
is performed by finding the face neighbors, i,e. neighbors which have one voxel face in 
5 common. As an alternative the surface of the vessel can be represented by surface triangles 
which can be generated by means of the so-called marching cubes algorithm which is known 
in the art The marching cubes algorithm is for instance described in the article: ^^Marching 
Cubes: A High Resolution 3D Surface Construction Algorithm'', by Lorensen and Cline, 
Computer Graphics, Vol. 21, No 4, July 1987. 
10 In a next step for each of the vessel surface points a gradient is calculated. 

This is perfoimed by means of standard techniques, which are among others described in the 
international patent application EP0Q/09S05. The resulting gradients for points 4, 5 and 6 are 
respectively indicated as v4» vS and v6. 

In Older to determine which of die three points should be included in the 
15 selection of points which are used for the calculation of the desired ordiogonal orientation of 
the cutting plane the mediod comprises the following characterizing steps. 

First for each point a vector is determined from the center of the tubular 
structure to said point For point 4, said vector is designated as vc4- For point 5, said vector is 
^g^gjiftfArf as vc5. For point 6, said vector is designated as vc6. 
20 Next for each surface point the angle is determined between each of said 

vectors and the gradient at said point. If said angle is equal to or smaller than a predetermined 
ceiling value said surface point is added to the above-mentioned selection of surface points, 
which are used to orient the cutting plane. A ceiling value can be set by a user or can be 
automatically set. Horn expmments a ceiling value between 50 and 70 degrees, for instance 
25 approximately 60 degrees, has proven to give reasonable results in the case of die current 
example relating to shape extraction for a blood vessel. 

Returning now to figure 1 we will perform some calculations in order to 
illustrate the method according to ±e invention. Starting with point 4 the angle between v4 
and vc4 is zero degrees. Let's assume that the ceiling value is chosen at 60 degrees. Point 4 is 
30 then added to die selection. This is rightly so because point 4 is a surface point on vessel 1, 
the vessel under analysis. However the angle betweo:! v5 and vc5 is more than 90 degrees. 
Point 5 will therefor under die current conditions be excluded from the selection. Again diis 
is conect, since point S is a surface point of vessel 2, the nei^boring vessel of the vessel 
investigated 
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The invention concerns a method of analyzing an object data set compiising 
points in a multi-dinxensional space, in which dacaset a tubular stmcnue occurs, said noethod 
compdsing tb& following steps: 

a) Choosing a starting position in or near the tubular structure; 
5 b) Deriving a cutting plane through tiie tubular structure at the starting 

c) Detennining a number of points forming part of the surface of the tubular structure in the 
vicinity of ^ starting position: 

d) Calculating a graxlient to the surface for each of said points. 
The n^thod furdier comprises the characterizing steps of; 

10 e) For each point deteimining a vector from the center of die mbular structure to said point; 

f) Determining the angle between said vector and the gradient at said point; 

g) Adding said point to a selection of points if said angle is equal to or smaller than a 
predetermined ceiling value; 

h) Using said selection of i>oints to calculate an orientation for the cutting plane such that 
IS the direction thereof is as parallel as possible to the longitudinal axis of the tubular 

structure at the starting position; 

i) Repeating steps a throu^ h for a new starting position along the tubular structure, if 
necessary. 

The invention also refers to a computw program to carry out the method 
20 according to the invention. 

Kg. 
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CLAIMS: 



1. A method of analyzing an object data set comprising points in a multi- 

dimensiooal space, in which dataset a tubular structuie occurs, said method comprising the 
following steps: 

a) Choosing a starting position in or near the tubular structuie; 

b) Deriving a cutting plane through the tubular structure at the starting position, 

c) Determining a number of points f oxming part of the surface of the tubular structure in the 
vicinity of the starting position; 

d) Calculating a gradient to the sur&ce for each of said points; 
characterized in that the method further comprises the steps of: 

e) For each point detemiining a vector from the center of the tubular stmcture to said point; 

f) Determining the angle between said vector and the gradient at said point; 

g) Adding said point to a selection of points if said angle is equal to or smaller than a 
predet^mined ceiling value; 

h) Using said selection of points to calculate an orientation for the cutting plane such that 
the direction thereof is as parallel as possible to tfie longitudinal axis of the tubular 
stmctuie at the stardng position; 

i) Reputing steps a through h for anew starting position along the tubular structuie, if 
necessary. 

2. Method according to claim l,fuxther comprising the st^s of: 
defming a sphere, which is at least partially intersected by the tubular structure; and 
pexfoiming steps e through g only for points lying inside die sphere. 

3. Method according to claim 1 or 2, wherein steps e throu^ g are pexfonned 
only for points lying at a piedetennined maximum distance from the cutting plane» 

4. A computer program to carry out the method according one or more of the 
preceding claims. 
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It is noted that the angle betwem v6 and vc6 appears to be smaller than 60 
degrees* Point 6 would unintendedly be added to the selection, since point 6 is a surface point 
of vessel 2. This can be pRV^ted by adding another criterion relating to the selecti<m of 
vertices. A useful second CTltenon is that only surface points lying inside sphere 3 should be 
5 included in the selection. From figure 1 it is clear that point 6 does not meet this criterion and 
would therefor not be included in the selection after all. 

It is noted that a useful third criterion would be to only perform steps e 
through g for points lying at a predetermined maximum distance from the cutting plane. 
When de object data set comprises surface triangles, a ceiling value for this distance can e.g. 
10 be denved form the length of the triangle edges (see also the above mentioned article by 
Bruijns). 

The order in which the three criteria according to tiie method accotding to the 
invention are implemented may be varied. 

The thus created selection of suxf ace points is now used to orient the cutting 

1 5 plane such tbat the directicm thereof is as parallel as possible to the longitudinal axis of the 
mbular structure at the starting position. An example of this calculation is described in the 
article mentioned above, which document is entirely incorporated herein by reference. 
The above steps of tiie method according to die invention can be repeated for each new ; ) 

Starting position along the tubular structure, if necessary. 

20 Now the method of the invention is explained, a skilled person will be able to 

translate the steps of the method into a computer program to carry out the method. 

Summarizing the invention refers to a method for analyzing an object data set 
in a mulii-dimensional space comprising a tubular structure. This method is especially useful 
in the field of medical diagnostics and treatment, where the object notably is a patient to be 

25 examined. A]l kinds of tubular structures can be analyzed, such as blood vessels, but also 
brain cell or neuron cell structures. In analyzing blood vessels the method according to the 
invffltion allows for an accurate detemiination of the size and diameter of a blood vessel 
which is essential for the ri^t diagnosis of e.g. stenosis or an aneurysm and the save 

treatment thereof. 

30 The invention is of course not limited to the described or shown embodiment, 

but generally extends to any embodiment, which falls within the scope of the appended 
claims as seen in li^t of the foregoing description and drawings. 



Ontvangst tijd 21. dec. 12:30 



Afdruk tijd 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the appHcant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
JzfpADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please dp not report these problems to 
the IFW Image Problem Mailbox. 



I: 



TH\SPAGEBlANK(us"0! 



